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ABSTRACT 


Enstatite (EngoFs^Wooi) was synthesized from a hypersthene normative basaltic melt 
doped at the same time with La, Ce, Nd, Sm, Eu, Dy, Er, Yb and Lu. The rare earth element 
concentrations were measured in both the basaltic glass and the enstatite. Rare earth element 
concentrations in the glass were determined by electron microprobe analysis with uncertainties 
less than two percent relative. Rare earth element concentrations in enstatite were determined by 
secondary ion mass spectrometry with uncertainties less than five percent relative. The resulting 
rare earth element partition signature for enstatite is similar to previous calculated and composite 
low-Ca pigeonite signatures, but is better defined and differs in several details. The partition 
coefficients are consistent with crystal structural constraints. 
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INTRODUCTION 


As part of on-going research into the petrogenesis of basaltic achondrite meteorites, we 
have recently determined mineral/melt partition coefficients for nine rare earth elements (REE) in 
enstatite using the percent-level doping technique. Hypersthene normative basaltic glass was 
concurrently doped with La, Ce, Sm, Nd, Eu, Dy, Er, Yb, and Lu. Isothermal growth 
experiments were conducted at 1210°C under oxygen fiigacity conditions equivalent to the iron- 
wustite buffer. These conditions produce enstatite crystals (EngoFsisWooi) with dimensions of 
approximately 75 x 150 x >1500 pm. The REE abundances were measured in the pyroxene with 
a Cameca 4f ion microprobe. REE abundances in the basaltic glass were measured by electron 
microprobe analysis using the Cameca Camebax and Cameca SX-100 microprobes at the Johnson 
Space Center. These results provide a REE partition coefficient signature for enstatite. These 
partition coefficients may be more appropriate for natural enstatite crystals than composite 
partition coefficients derived for low-Ca pigeonite (e.g., McKay et al., 1986, 1991), because in 
the present experiments (1) the pyroxene is enstatite rather than pigeonite, and (2) all nine REE 
were present simultaneously in each charge rather than as individual REE in different charges as in 
the earlier experiments, and were therefore subject to identical experimental and analytical 
conditions. This approach provides more assurance of a self-consistent REE partition coefficient 
signature than calculated or composite signatures derived from several singly doped experiments, 
because experimental and analytical conditions are sometimes difficult to reproduce exactly, and 
in singly doped experiments the dopant elements interact in a more simple chemical system. 
Therefore we present enstatite partition coefficients for nine of the fifteen rare earth elements. 
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EXPERIMENTAL PROCEDURE 


Five pure oxides; Si02, AI2O3, MgO, Fe203, and CaO were mixed in appropriate 
proportions to produce 10 g of a basaltic glass with a composition similar to the basaltic 
achondrite meteorites, eucrites, though somewhat more hypersthene normative. The oxide 
mixture was ground for half an hour under acetone in a motorized mortar and pestle, then melted 
at 1350 °C in a large platinum crucible in air for 24 hours. The melt was quenched in de-ionized 
water and then crushed and ground in the mortar and pestle under acetone for about fifteen 
minutes to produce a fine glass powder. Cr 2 03 , TiC> 2 , La, Ce, Nd, Sm, Eu, Dy, Er, Yb oxides and 
Lu(NC> 3)3 were weighed out and added to the appropriate amounts of glass powder and then 
thoroughly mixed under acetone with a mortar and pestle. Total minor and trace element 
concentrations were held to less than 4 wt%. 

Mineral growth experiments were initiated by suspending 125 mg pressed glass-powder 
pellets on 3 mm diameter platinum wire loops in a flowing CO/CO 2 gas-mixing furnace at 1300 
°C and iron-wiistite redox conditions for 48 hours to equilibrate the samples with the fugacity of 
the experiments and to equilibrate iron exchange with the platinum loop. Approximately 5 wt% 
additional Fe 2 C >3 was added to the original mixture to compensate for iron loss to Pt in all stages 
of the experiments. Because of the much higher diffusion rates at 1300°C than at the 
crystallization temperature of 1210°C, most iron loss occurred during this high temperature 
equilibration. After equilibration, the samples were air quenched by rapid removal from the 
furnace. Specific experiments were then initiated by placing a sample back into the furnace at 
iron-wiistite fugacity conditions, superheating it to 1260°C (~50 degrees above the liquidus). 
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holding it there for ~4 minutes, and then dropping the furnace temperature rapidly to 1210°C (2 ± 
1 degrees below the liquidus temperature). The samples were then held at constant temperature 
for 120 hours, and air-quenched by pulling the sample support assembly out of the furnace. This 
procedure produces samples having less than five percent total crystallization and minimizes 
nucleation, whereby the few crystals that grow are large enough to analyze without difficulty (75 
x 150 by >1000 pm). Nucleation appears to be internal to the sample charge, as the crystals 
predominantly float freely in the basaltic glass droplet and do not appear to nucleate on the inside 
surface of the glass sphere or on the wire loop (Lofgren, 1980). 

Experiment temperatures are monitored with a thermocouple mounted in the furnace 
muffle tube approximately three millimeters from the sample. The Pt and Pt-Rh thermocouple is 
calibrated against the melting point of gold (1064 °C). Temperatures are believed to be accurate 
within ± 5 degrees, with relative precision of about ± 1 degree. The oxygen fiigacity is monitored 
using a solid-state oxygen sensor in a separate reference furnace (Jurewicz et aL, 1993) calibrated 
using a commercially mixed 1 : 1 CO:C02 reference gas. The accuracy for the oxygen fiigacity 
measurements is better than ± 0. 1 log units, and the experimental oxygen fiigacity was equivalent 
to iron-wiistite minus 0.2 log units ±0.1 log units (log/ O 2 = -1 1.8 at 1210 °C). 

ANALYTICAL PROCEDURES 

Major and minor element compositions of the pyroxene and glass were determined with a 
five spectrometer Cameca SX-100 electron microprobe. Quantitative analyses utilized natural 
mineral standards, 15 kV, 40 nA, and the Cameca PAP correction routine (Pouchou and Pichior, 
1991). One-sigma counting uncertainties for the major elements are less than one percent of the 
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amount present. 


Rare earth element concentrations of the enstatite crystals were determined using the 
University of New Mexico-Sandia National Laboratories, Advanced Materials Laboratory, 
Cameca 4f ion microprobe. The 20 nA primary O minus beam was focused to a 20-25 pm beam 
diameter. The dopant concentrations in these synthetic pyroxenes yield secondary ion count rates 
high enough to provide theoretical standard deviations of less than one percent. However, due to 
uncertainties in the standard compositions, the analytical uncertainties are conservatively 
estimated at less than 5 percent. The nine REE used in these experiments were chosen to avoid 
both x-ray spectra interferences and isobar interferences during secondary ion mass spectrometry. 
Ion probe analyses utilized the following isotopes: 139 La + , 140 Ce + , 14< Nd + , 147 SnT, l51 Eu + , l53 Eu + , 
163 Dy + , 167 Er + , 174 Yb + , and l75 Lu + . Enstatite REE analyses were calibrated against two natural 
pyroxene standards. The basaltic glass compositions were determined by electron microprobe 
analysis using a Cameca Camebax and a Cameca SX-100 microprobe, and synthetic glasses 
containing rare earth elements as standards (Drake and Weill, 1972). The counting uncertainties 
for the glass analyses are less than two percent. 

RESULTS AND DISCUSSION 

The synthesized pyroxene crystals have homogeneous major element compositions 
consistent with enstatite (Table 1). Longhi and Boudreau (1980) and Longhi (1987) showed that 
orthoenstatite was stable to “at least as low as 121 1 °C and presumably much lower. . .” 
Interestingly, the enstatite crystals of this study are sector zoned with respect to Al, Cr, and to a 
lesser extent the REE. Abundances of heavy REE show positive correlations with Al abundance 
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from one sector to another (Fig. 1,2), though the light REE show very little correlation with Al, 
as indicated by the low value of R 1 value for La as compared with Lu (Fig. 2). The heavy REE 
have more mild negative correlations with Mg and Fe. Detailed discussion regarding the 
development of minor and trace element sector zoning is in a manuscript in preparation by 
Schwandt and McKay. 

We define the partition coefficients ( D ) presented here according to the terminology of 
Beattie et al. (1993), where mBenVmelt £) M = maml C]J mcit Cu. M represents the cation of interest and 
C represents the weight fraction of the cation in each respective phase. The D ' s determined in 
each sector for the measured REEs are reported in Table 2. Figure 3 represents a plot of the 
partition coefficients against the REE’s in order of increasing atomic number. 

The resulting REE partition coefficient signature (Fig. 3) is similar to the theoretical 
signature as calculated from pyroxene wollastonite content with the equation of McKay et al. 
(1986) and to the composite low-Ca pigeonite signature (McKay et al. 1991). However, because 
more REE were measured, the new partition coefficient pattern is better defined than the older 
ones, especially for the HREE. The new pattern also differs in important details from the older 
ones. Compared with the composite pattern from McKay et al. (1991), the new pattern has a 
steeper slope for the LREE, and a higher Dsm/D N( i ratio. This may be a result of the difficulties of 
measuring very low partition coefficient values using the electron microprobe (McKay, 1986). 
Compared with the pattern calculated using the equation of McKay et al. (1986), the new pattern 
is steeper overall, and much smoother. This result is not surprising because the calculated pattern 
in Fig. 3 represents an extrapolation to significantly lower Wo content than experimental 
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pyroxenes used by McKay et al. (1986) to parameterize ^Dree as a function of Wo. 

The new partition coefficient pattern also shows a significant negative Eu anomaly, 
indicating that the partition coefficient for Eu +2 is less than that for Eu +3 . Of course the size of this 
Eu anomaly will vary with oxygen fugacity, with the current experiments having been run at 
fugacities close to IW. Thus, at low /oj, orthopyroxene will have a negative Eu anomaly even if it 
crystallizes from a melt with chondritic Eu/Sm. 

The observed REE partition coefficient signatures are easily rationalized based on crystal 
structural constraints. The size of the pyroxene M2 site is more flexible than the Ml site, though 
cations larger than Ca are not easily incorporated into the structure (Cameron and Papike 1981). 
Therefore, as the light REE are significantly larger in ionic size than the heavy REE it is logical 
that the light REE are depleted in enstatite relative to the heavy REE. The presence of a negative 
Eu anomaly indicates that the larger ionic radius of Eu +2 (1.25A in eight-fold coordination, vs. 
1.07 A for Eu +3 , Shannon, 1976) must make substitution more difficult than for Eu +3 even though 
the latter requires a charge-balancing coupled substitution. 

Based on pyroxene compositions from lunar and meteoritic basaltic rocks, Papike (1980) 
identified a few primary substitution mechanisms for minor elements in pyroxene: 1) ^Al 3 * + 
1V A1 3+ for VI Mg 2+ + 2) VI Cr 3+ + IV A1 3+ for ^Mg 2 " + w Si*. 3) ^Ti 4 * + 2 IV A1 3+ for ^Mg 2 * + 

IV Si 4+ . Similarly, relative to the sector zoning in our synthetic pyroxenes, Mg and Si show strong 
inverse correlations and Cr a strong normal correlation with Al concentration. This suggests the 
same coupled substitution mechanisms are important. Incorporation of the REE into the enstatite 
structure must also be occurring via a similar coupled substitution ( VI REE 3+ + IV A1 3+ for VI M2 2+ + 
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^Si 4 *) because REE concentrations in the enstatite correlate positively with A1 concentration (at 
least for the HREE, Fig. 2) and correlate negatively with the typical M2 cations. The REE 
correlations are strongest with A1 and only moderate relative to the M2 cations individually. 
Therefore this REE partitioning behavior suggests that the incorporation of impurity cations, such 
as the REE, in the enstatite crystal structure is at least partly controlled by the substitution of A1 
for Si in the tetrahedral sites. Moreover, these results suggest there is an additional charge 
compensation involving the larger impurity cations in the octahedral M2 site. The mechanism 
behind this observation and investigation of the mechanisms leading to the development of sector 
zoning of minor and trace elements is explored in depth in a manuscript in preparation by 
Schwandt and McKay. 

CONCLUSION 

We have determined nine REE partition coefficients and the resulting REE partitioning 
signature for enstatite synthesized from a hypersthene normative basaltic melt. We believe these 
partition coefficients provide the most complete and reliable experimentally determined partition 
coefficient pattern currently available for orthopyroxene. These partition coefficients are useful 
for research on the petrogenesis of basaltic achondrite meteorites and terrestrial mantle materials. 
In addition, the results also provide clues about mechanisms resulting in the incorporation of 
impurity cations into the pyroxene crystal structure. 
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Table 1. Major element compositions determined by electron 
microprobe analysis. Trace element values in glass determined 
by electron microprobe analysis. Trace element values in 
enstatite obtained by SIMS. Column numbers match SIMS 
analysis spots shown in Figure 1 . 



Glass 

Enstatite 




wt% 


1 

2 

3 

4 

5 

Si0 2 

53.31 

52.96 

53.26 

53.79 

53.33 

53.44 

A1 2 0 3 

15.07 

3.29 

3.42 

2.52 

3.30 

3.01 

MgO 

8.06 

28.37 

28.49 

28.66 

28.53 

28.35 

FeO 

12.50 

12.64 

12.84 

12.99 

12.77 

12.90 

CaO 

6.93 

0.63 

0.64 

0.63 

0.64 

0.64 

Ti0 2 

0.22 

0.05 

0.05 

0.05 

0.06 

0.06 

Cr 2 0 3 

0.13 

0.92 

0.98 

0.78 

0.86 

0.84 

ppm: 







La 

5012 

1.99 

1.62 

1.43 

2.82 

2.65 

Ce 

6118 

6.72 

5.74 

4.74 

7.5 

8.49 

Nd 

4322 

18.5 

16.18 

12.56 

17.68 

21.39 

Sm 

4831 

59.6 

54.0 

40.0 

57.5 

65.6 

Eu 

4738 

21.9 

21.3 

15.8 

21.7 

19.8 

Dy 

3236 

171 

175 

121 

173 

178 

Er 

3297 

305 

320 

225 

323 

317 

Yb 

3664 

381 

412 

294 

419 

394 

Lu 

4641 

478 

533 

384 

532 

503 

Mole Fraction: 






fe' 

0.46 

0.19 

0.20 

0.20 

0.20 

0.20 

Wo 


0.012 

0.012 

0.012 

0.012 

0.013 

En 


0.79 

0.79 

0.79 

0.79 

0.80 


fe' = FeO/(FeO+MgO); Wo = CaO/(CaO+FeO+MgO); 
En = MgO/(CaO+FeO+MgO); all quantities molar. 



Table 2. Rare earth element partition coeffi- 
cients. Column numbers match SIMS analy- 
sis spots shown in Figure 1 . Uncertainties 


less than 5%. 


1 


La 

0.0004 

Ce 

0.0011 

Nd 

0.0043 

Sm 

0.0123 

Eu 

0.0046 

Dy 

0.052 

Er 

0.092 

Yb 

0.104 

Lu 

0.103 


2 

3 

0.0003 

0.0003 

0.0009 

0.0008 

0.0037 

0.0029 

0.0112 

0.0083 

0.0045 

0.0033 

0.054 

0.037 

0.097 

0.069 

0.112 

0.080 

0.114 

0.082 


4 

5 

0.0006 

0.0005 

0.0012 

0.0014 

0.0041 

0.0049 

0.0119 

0.0136 

0.0046 

0.0042 

0.053 

0.054 

0.098 

0.096 

0.114 

0.107 

0.114 

0.108 
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Figure Captions: 

1 . Digital map of A1 concentration. Bright indicates high A1 concentration and dark indicates 
low A1 concentration. Enstatite crystal with sector zoning of A1 surrounded by basaltic glass. 
Circles outline ion microprobe analysis locations. Analyses 1 and 5 are in { 100 } sectors, 2 
and 4 are in {010} sectors and 3 is in a portion of {001 } sector overlapped by a low A1 zone 
that is the result of growth twinning. 

2. Plot of AI 2 O 3 wt % versus the REE concentration (ppm) showing the correlation of the REE 
with Al. Correlation Coefficients (R 2 ) are given for La and Lu, and indicate a declining 
correlation with increasing cationic size, i.e. lighter REE. 

3. Partition coefficient versus atomic number plot reveals the rare earth element partitioning 
signature for enstatite. The field represents the variation of the partition coefficients as a 
function of the sector zoning. The experimental signature is compared with a pattern 
calculated from the variation of partition coefficients with the wollastonite content of synthetic 
shergottite pyroxenes (McKay et al., 1986) and a composite low-Ca pigeonite pattern from 
several experiments doped with single REE (McKay et al., 1991). 
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